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Abstract: In the framework of the CMSSM we study the gravitino as the hghtest super- 
symmetric particle and the dominant component of cold dark matter in the Universe. We 
include both a thermal contribution to its relic abundance from scatterings in the plasma 
and a non-thermal one from neutralino or stau decays after freeze-out. In general both 
contributions can be important, although in different regions of the parameter space. We 
further include constraints from BBN on electromagnetic and hadronic showers, from the 
CMB blackbody spectrum and from collider and non-collider SUSY searches. The region 
where the neutralino is the next-to-lightest superpartner is severely constrained by a con- 
servative bound from excessive electromagnetic showers and probably basically excluded 
by the bound from hadronic showers, while the stau case remains mostly allowed. In both 
regions the constraint from CMB is often important or even dominant. In the stau case, for 
the assumed reasonable ranges of soft SUSY breaking parameters, we find regions where the 
gravitino abundance is in agreement with the range inferred from CMB studies, provided 
that, in many cases, a reheating temperature Tr is large, ~ 10^ GeV. On the other 
side, we find an upper bound Tr < 5 x 10^ GeV. Less conservative bounds from BBN or an 
improvement in measuring the CMB spectrum would provide a dramatic squeeze on the 
whole scenario, in particular it would strongly disfavor the largest values of Tr ~ 10^ GeV. 
The regions favored by the gravitino dark matter scenario are very different from standard 
regions corresponding to the neutralino dark matter, and will be partly probed at the LHC. 
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1. Introduction 

Low-energy super symmetry (SUSY) provides perhaps the most attractive candidates for 
cold dark matter (CDM) in the Universe. This is because in the SUSY spectrum several 
new massive particles appear, some of which carry neither electric nor color charges. The 
lightest among them (the lightest SUSY partner, or the LSP) can then be neutral and 
either absolutely stable by virtue of some discrete symmetry, like i?-parity, or very long- 
lived, much longer than the age of the Universe, and thus effectively stable. A particularly 
well-known and attractive example of such a weakly-interacting massive particle (WIMP) 
is the lightest neutralino. In recent years, however, there has been also a renewed interest 
in the two alternative well-motivated SUSY candidates for WIMPs and CDM, namely the 
gravitino and the axino. 

The spin-3/2 gravitino acquires its mass from spontaneous breaking of local SUSY, or 
supergravity. Since its interactions with ordinary matter are typically strongly suppressed 
by an inverse square of the (reduced) Planck mass, cosmological constraints become an is- 
sue [1]. Early on it was thought that, with a primordial population of gravitinos decoupling 
very early, if stable, they had to very light, below some < IkeV [2], in order not to over- 
close the Universe, or otherwise very heavy, > lOTeV [3], so that they could decay before 
the period of BB nucleosynthesis (BBN). With inflation these bounds disappear [4, 5] but 
other problems emerge when gravitinos are re-generated after reheating. If the gravitino 
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is not the LSP, it decays late (~ 10*^ sec) into the LSP (say the neutralino) and an ener- 
getic photon which can distort the abundances of hght elements produced during BBN, for 
which there is a good agreement of calculations with direct observations and with CMB 
determinations. Since the number density of gravitinos is directly proportional to the re- 
heating temperature Tr, this leads to an upper bound of Tr < 10^~^ GeV [4, 5, 6, 7, 8, 9] 
(for recent updates see, e.g., [10, 11]). On the other hand, when the gravitino is the LSP 
and stable, ordinary sparticles can decay into it and an energetic photon. A combination 
of this and the overclosure argument (Og/i^ < 1) leads in this case to an upper bound 
rR<109GeV [6, 12]. 

Because of many similarities, for comparison we comment here on axinos as DM. 
The axino is a fermionic superpartner of the axion in supersymmetric models with the 
Peccei— Quinn (PQ) mechanism implemented for solving the strong CP problem [13]. Axino 
interactions with ordinary matter are suppressed by 1/ fa, where fa ~ 10^^ GeV is the PQ 
scale. In many models the axino mass is not directly determined by a SUSY breaking 
scale MsusYj in contrast to the neutralino and the gravitino, and therefore the axino can 
naturally be the LSP. Without inflation the axino has to be light (< 1 keV) and thus warm 
DM [14, 15, 16, 17, 18]. Otherwise it can naturally be a cold DM [19, 17, 20], so long as 
Tr < 10^~^ GeV. Constraints from Big Bang Nucleosynthesis (BBN) on axino CDM are 
relatively weak since NLSP decays to them typically take place before BBN. See [21] for 
an improved treatment of thermal production of axinos at Tr > 10^ GeV. 

For both the gravitinos and the axinos, despite their typical interaction strengths being 
so much weaker than electroweak, their relic abundance can still be of the favored value of 
~ 0.1. This is because they can be efficiently produced in a class of thermal production (TP) 
processes involving scatterings and decays of particles in the primordial plasma, depending 
on Tr. Alternatively, in a non-thermal production (NTP) class of processes, the next -to- 
lightest supersymmetric particle (NLSP) first freezes out and next decays to the axino or 
gravitino. These mechanisms are supposed to re-generate the relics, after their primordial 
population has been diluted by a proceeding period of inflation. 

In TP, gravitino (or, alternatively, axino) production proceeds predominantly through 
ten classes of processes involving gluinos [5, 12]. In four of them, a logarithmic singularity 
appears due to a t-channel exchange of a massless gluon which can be regularized by 
introducing a thermal gluon mass [5, 12]. A full result for the singular part was obtained 
in [22]. In [23] a resumed gluon propagator was used to obtain the finite part of the 
production rate, and an updated expression for the relic abundance fl^h^ of gravitinos 
generated via TP, valid at high Tr, was given 

n-Z^h^ ~ 2 ( ^ ( 100 Gey \ (rn^\ ' 

where mg{^) above is the running gluino mass. In [22, 23] it was argued that, for natural 
ranges of the gluino and the gravitino masses, one can have Vl~^h^ ~ 0.1 at Tr as high as 
-[^q9-io QgY_ Such high values of Tr are essential for thermal leptogenesis [24, 25], with a 
lower limit of Tr > 2 x 10^ GeV [26]. 
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The issue of gravitino relics generated in NTP processes and associated constraints 
was recently re-examined in detail in [27, 28, 29, 30] and in [31]. Since all the NLSPs 
decay into gravitinos, in this case 

nfh' = ^^^^NLSP^^ (1.2) 

where Q^'^^h^ is the NTP contribution to the gravitino relic abundance and ^nlsp^^ 
would have been the relic abundance of the NLSP if it had remained stable. Note that 

grows with mp^. 

In grand-unified SUSY frameworks, like the popular Constrained Minimal Supersym- 
metric Standard Model (CMSSM) [32], which encompasses a class of unified models where 
at the GUT scale gaugino soft masses unify to mi/2 and scalar ones unify to mo, one 
often finds that the relic abundance of the neutralino (or the stau) is actually significantly 
larger than 0.1. One therefore in general cannot neglect the contribution from the NTP 
mechanism, unless mg is rather small. In this case, however, at high Tr ~ 10^ GeV, TP 
contribution is likely to play a role, since oc T^/uIq. Clearly, in general both TP 

and NTP must be simultaneously considered. 

In [30] NTP of gravitinos was considered in an effective low-energy SUSY scenario. 
The relic abundance of gravitinos from NTP via neutralino, stau (which had already been 
examined in gauge— mediated SUSY breaking schemes in [33]) and sneutrino NLSP decays 
was however only crudely approximated and that from TP was not included at all. A 

weak constraint VlJi^^h'^ < 0.1 (rather than Q^'^'^h'^ ~ 0.1) was assumed. On the other 

G G 

hand, constraints from BBN were treated with much care. Typically, lifetimes for NLSP 
decays into gravitinos are ~ 10^ sec in which case constraints from electromagnetic fluxes 
are particularly important. Nevertheless, hadronic showers, which in the past were thought 
to be important only for lifetimes < 10^ sec, must also be included in considering particle 
decays in late times, since they provide additional strong constraints. In the case of light 
gravitinos in gaugc-nicdiatcd SUSY breaking schemes this constraint was applied in [34] 
and in the case of CDM gravitinos in [29, 30]. Furthermore, in [29, 30] substantial con- 
straints on the SUSY parameter space were derived from a requirement of not distorting 
the CMB blackbody spectrum by energetic photons [1, 5]. 

NTP of gravitinos in the framework of the CMSSM was examined in [31]. Constraints 
from electromagnetic showers were applied, but not from hadronic ones. Nor was the 
constraint from CMB applied. Gravitino abundance from NTP was computed much more 
accurately than in [27, 28, 29, 30]. On the other hand, similarly to [27, 28, 29, 30], only cos- 
mologically allowed regions VL^^^h? < 0.1 were delineated but not cosmologically favored 
ones of Q^^h? ~ 0.1. For the assumed ranges of parameters some regions were found 

G 

where Q^^h? was not excessively large, but actually too low. (In [31] it was also noted 
that any possible stau NLSP asymmetry would be washed away by stau pair-annihilation 
into tau pairs.) As we will show later, unlike the authors of [31], in the CMSSM at large 
values of mi/2 (beyond those considered in [31]) we have found cosmologically favored 
regions of Q^^^h? ~ 0.1. 
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A question arises for which values of SUSY parameters, as weh as Tr, the combination 
of gravitino yields from both TP and NTP gives O^/i^ ~ 0.1 in unified SUSY schemes. In 
this paper we investigate this issue within the CMSSM, which is a model of much interest. 
We assume no specific underlying supergravity model and treat as a free parameter. 
In the CMSSM soft masses are assumed to be generated via a gravity-mediated SUSY 
breaking mechanism, in which case mi/2; "j-o can be in the GeV to TeV range, 

and we need to ensure that tlic gravitino is the LSP. We compute the relic abundance 
of gravitinos with high accuracy which matches present observational precision of CDM 
abundance determinations. In evaluating O.J^h'^ we follow [23], while OM^^^h'^ is determined 
by the yield of the NLSP which we compute numerically, following our own calculation 
(without partial wave expansion), as described below. Wc apply constraints from both the 
electromagnetic and from the hadronic fluxes and from CMB spectrum, as well as the usual 
constraints from collider and non-collider SUSY searches. We concentrate on the largest 
Tr, ~ 10^ GeV but also consider the region of low Tr where NTP dominates. 

Before proceeding, we should note that there are other possible gravitino production 
mechanisms, e.g. via inflaton decay or during preheating [35, 36], but they are much more 
model dependent and not necessarily efficient [37]. Alternatively, gravitinos may be pro- 
duced from decays of moduli flelds [38]. In this paper, we do not include these effects. 

Wc further assume i?-parity conservation, both for simplicity and because otherwise 
it is hard to understand why weak universality works so well. However, it is worth re- 
membering that in the case of such super-weakly interacting relics as the axino or the 
gravitino, i2-parity is not really mandatory, unlike in the case of the neutralino WIMP. 
Indeed, the suppression provided by the PQ or Planck scale is often sufficient to ensure 
effective stability of such relics on cosmological time scales even when i?-parity breaking 
terms are close to their present upper bounds. Indeed, in the case of the axino CDM, a 
tiny amount of its decay products into e"'"e~ pairs has been proposed as an interesting way 
of explaining an apparent INTEGRAL anomaly [39]. 

In the following, we will first summarize our procedures for computing fl^h^ via both 
TP and NTP. Then wc will list NLSP decay modes into gravitinos, and discuss constraints 
on the CMSSM parameter space, in particular those from BEN and CMB. Finally, wc will 
discuss implications of our results for thermal leptogenesis and for SUSY searches at the 
LHC. 

2. Gravitino Relic Abundance 

The present relic abundance of any stable, massive relics (x, f, d, G, . . .) produced either 
thermally or non-thermally is related to their yield^ as 



^We define the yield as Y = n/s, where s = (27r^/45)(;*sT^ is the entropy density, following a common 
convention of [40] which is also used in [11]. Another definition, used, e.g., in [30, 10], is Y' = n/n^ where 
n-y is the number density of photons in the CMB, n-^ = 2nrad = 2^(3)T^/7r'^. At late times t > 10® sec, 
typical for NLSP decays to gravitinos, and later, s ~ 7.04n-y. 




(2.1) 
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where m is the mass of the rehc particle. 
2.1 TP 

The yield of massive relics generated through TP processes can be obtained by integrating 
the Boltzmann equation with both scatterings and decays of particles in the expanding 
plasma [40]. In the case of the gravitino LSP, dominant contributions come from 2-body 
processes involving gluinos [5, 12, 22, 23]. These are given by a dimension-5 part of the 
Lagrangian describing gravitino interactions with gauge bosons and gauginos.^ For the ten 
classes of scattering processes the cross section, at large energies, has the form 



where Mp = I/a/SttG^ = 2.4 x 10^^ GeV is the reduced Planck mass. 

In actual computations we solve the Boltzmann equation numerically by following the 
usual steps described, e.g., in [23, 17], and use the expression (44) of [23] for the sum of 
soft and hard contributions to the collision terms. 

We do not include gluino decays into gravitinos which, like in the case of axino LSP, 
would only become important at Tr ~ nig [12, 17]. This is because we concentrate on large 
Tr ~ 10^ GeV which are relevant for models of thermal leptogenesis. 



In computing the relic abundance of gravitinos generated through NTP processes we 
first compute their yield after freezeout. In the CMSSM in most cases the NLSP is 
either the (bino-like) neutralino or the lighter stau. In the case of the neutralino, we 
include exact cross sections for all the tree-level two-body neutralino processes of pair- 
annihilation [42, 43] and coannihilation with the charginos, next-to-lightest neutralinos [44] 
and sleptons [45]. This allows us to accurate compute the yield and Q-^h^ in the usual case 
when the lightest neutralino is the LSP. We further extend the above procedure to the 
case where the NLSP is the lightest stau fi (a lower mass eigenstate of tr and fx,). We 
include all slepton-slepton annihilation and slepton-neutralino coannihilation processes. 
In both cases we numerically solve the Boltzmann equation for the NLSP yield and use 
exact (co) annihilation cross sections which properly take into account resonance and new 
final-state threshold effects. The procedure has been described in detail in [46] and was 
recently applied to the case of axino LSP [47]. 

3. NLSP Decays into Gravitinos 

Once the NLSPs freeze out from thermal plasma at t ~ 10~^^ sec, their comoving number 
density remains basically constant until they start decaying into gravitinos at late times 
i ~ 10^ - 10^ sec (so long as mg > O.lMgusY which we assume here). Associated decay 
products generate energetic fluxes which are mostly electromagnetic (EM) but also hadronic 





(2.2) 



2.2 NTP 



^See, e.g., [41]. 
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(HAD). If too large, these will wreck havoc on the abundances of light elements. Limits 
on electromagnetic showers become important for t > 10"^ sec and come mainly either 
from excessive deuterium destruction via 7 + Z> — n + p (10^ sec < t < 10^ sec), or 
production via 7 + ^He D + . . . > 10^ sec). Decay rates and branching ratios into EM 
radiation generated by late-decaying particles have recently been re-analyzed by Cyburt, 
et al, in [10]. Updated bounds on EM fluxes have been obtained on the parameter (^x = 
iTT'Xnx/n.y, where X denotes the decaying particle,^ rij is a number density of background 
photons, as a function of X lifetime tx, by assuming that in X ^ G + . . . decays associated 
showers are mostly electromagnetic. This is indeed the case when X is either the neutralino 
or the stau [31, 29]. 

Limits from BBN on hadronic showers are stronger for r < 10^ sec but are often also 
important at later times [48, 49, 11, 50]. They come mainly from '^He overproduction via 
n + p ^ D '^He for tx < 10^ sec, and D overproduction via n + '^He ^ D at later 
times. Hadronic components produced in late X decays into gravitinos, while much less 
frequent, will still lead to important constraints on the parameter space, as mentioned 
above, and will play an important role in our analysis. Upper limits on hadronic radiation 
from X decays, but with somewhat stronger assumptions than in [10], have recently been 
re-evaluated by Kawasaki, et al, in [11]. 

More stringent bounds, by roughly a factor of ten, come from considering constraints 
from %i [51] and/or from '^He [52]. As discussed, e.g., in [29] at present both are probably 
still too poorly determined to be treated as robust. For this reason, and in order to remain 
conservative, we do not use the constraint from ^Li, unlike in [31]. Nor, like in [31], do we 
apply the constraint from ^He. 

In order to apply bounds from BBN light element abundances on EM/HAD showers 
produced in association with gravitinos, we need to evaluate the relative energy (i = 
em, had), as defined below, which is released in NLSP decays into EM/HAD radiation. 
First, for each NLSP decay channel we need to know the energy ef transferred to EM/HAD 
fluxes. In decays X G + R + . . ., where R collectively stands for all the particles 
generating cither EM or HAD radiation, the total energy per NLSP decay carried by R 
will be a fraction of mx. This is because, at late times of relevance to G production, the 
NLSPs decay basically at rest. In 2-body decays 

where now ttir stands for mass of R. Unless mg is not much less than mx, then, for 
negligible mji, the usual approximation ~ W'x/2 works well. In the case of 3- and 
more-body final states E^^^ as a fraction of mx assumes a range of values. 

We also need to compute the NLSP lifetime tx and branching fractions Bf {i = 
em, had) into EM/HAD showers. All the above quantities depend on the NLSP and (with 
the exception of the yield) on its decay modes and the gravitino mass. For the cases of 
interest {x and fi) these have been recently evaluated in detail in [29, 30] (see also [31]) 
and below we follow their discussion. 



' Hereafter we will mean X = x,ti for brevity. 
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For the neutralino NLSP the dominant decay mode is % — > G7 for which the decay 
rate is [31, 30] 




where x = ^iiB + Ni2W^ + A^is^Tj + Ni^H^. In the CMSSM the neutrahno is a nearly 
pure bino, thus x — B. The decay x produces mostly EM energy. Thus 

B^m ^ 1, (3-3) 

2 2 
mr, — m ~ 

e^m = ^ (3-4) 

and the energy ^em released into electromagnetic showers is in this case simply given by 

CX fX TJX ("^ 

If kinematically allowed, the neutralino can also decay via x ~^ GZ, Gh, GH, GA for 
which the decay rates are given in [31, 30]. These processes contribute to hadronic fluxes 
because of large hadronic branching ratios of the Z and the Higgs bosons {B^^^ ~ 0.7, 
Bhad — O-^)- ^^^^ energy ^^^^ released into hadronic showers is 

^Ld^{Y.h^LA)Y'' (3-6) 
where the sum goes over all hadronic decay modes. 



V- X ox ^ e|r(x ^ GZ)B^^, + E elTjx ^ Gh)Bl, + e^.Fjx ^ Gqq) 

^ r (x - G7) + r (x - Gzj + r (x -> Gh) 



(3.7) 



where 

ml - m% + ml 

4 ~ ^' ^""^ k = Z,h, H, A, (3.8) 

X 

2 

egg ~ 3("ix -"^g)- (3-9) 

Below the kinematic threshold for neutralino decays into G and the Z/Higgs boson, 

one needs to include 3-body decays with the off-shell photon or Z decaying into quarks 
for which B^^^{x G'^* jZ* Gqq) ~ 10^^ [30]. This provides a lower bound on B^^^. 
At larger m^^ Higgs boson final states become open and we include neutralino decays to 
them as well. 

The dominant decay mode of the stau fi is fi Gt for which (neglecting the tau- 
lepton mass) the decay width is [31, 30] 

1 m%_ 



(~ \ m- I lib- \ 

fi ^ Gt] = 1 - ^ . (3.10) 
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In [28, 29, 30] it was argued that decays of staus contribute basically only to EM 
showers, despite the fact that a sizable fraction of tau-leptons decay into light mesons, like 
pions and kaons. These decay electromagnetically much faster than the typical time scale 
of hadronic interactions, mainly because at such late times there are very few nucleons left 
to interact with [28]. Thus 

Blln ^ 1, (3.11) 
2 2 

1 "^fi ~ 

2 2mf, ' ^^-^^^ 

where the additional factor of 1/2 appears because about half of the energy carried by the 
tau-lepton is transmitted to final state neutrinos. The energy released into electro- 
magnetic showers is in this case 

Qn^^7mBll^Y^K (3.13) 

As shown in [29], for stau NLSP, the leading contribution to hadronic showers come 
from 3-body decays fi GtZ^Gvt-W , or from 4-body decays fi — > Gt^*/Z* Grqq. 
The corresponding energy ^^^^ is 

Cd^(E4Xad)^'^ (3-14) 

where the sum goes over all hadronic decay modes, 

V- ^, e^Tjn ^ GrZ)Bl, + e%T{n ^ Gu^W)BZ, + e%T{n ^ Grqq) 

^ r (fi ^ Gr) 

and 

- - 41 ~ ^(^n - (3.16) 

One typically finds [30] ~ 10~^ — 10~^ when 3-body decays are allowed and ~ 10~® 
from from 4-body decays otherwise, thus providing a lower limit on the quantity. Given 
such a large variation in S^^^, the choice (3.16) is probably as good as any other. 



4. Constraints 

• The relic abundance We will be mostly interested in the cases where the sum 
VIqU^ = VL^h'^ + Q^^h"^ satisfies the 2a range for non-baryonic CDM 

0.094 < ^Qh^ < 0.129, (4.1) 

which follows from combining WMAP results [53] with other recent measurements of the 
CMB. Larger values are excluded. Lower values are allowed but disfavored. We will 
also delineate regions where alone satisfies the range (4.1). These regions will be 

q 

cosmologically favored for Tr <C 10 GeV when TP can be neglected. 
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• Electromagnetic and hadronic showers and the BBN As stated above, bounds 
on electromagnetic fluxes have recently been re-evaluated and significantly improved in [10], 
while on hadronic ones in [11]. (A clear summary of the leading constraints can be found 
in [29].) As mentioned above, each analysis uses somewhat different assumptions and input 
parameters. 

In constraining EM showers, Cyburt, et al, [10] imposed the following observational 
bounds on light element abundances 

1.3 X 10"^ < D/H < 5.3 X 10~^ (4.2) 
0.227 <Yp< 0.249, (4.3) 
9.0 X 10~^^ < (^Li/H) < 2.8 x 10"^°. (4.4) 

In the ^He abundance (Yp) the error bars were taken at 2 a. 

In deriving constraints from EM showers, we impose < (Cx/2) (n^/s) where (x 
is defined above and in [10] and the factor of n-y/s ~ 1/7.04 is due to a difference in 
our definitions of the yield. The largest allowed values of (^x, as a function of tx, for 
the assumed abundances of D/H + Yp + "^Li can be read out of fig. 7 of [10]. We fix 
the baryon-to-entropy ratio r] at 6.0 x 10~^^, which is consistent with the WMAP result 
= Q.ltoi X 10-^° [54]. 

In constraining hadronic fluxes from NLSP decays we compare with Kawasaki, et 
al., [11], who used the following input 

2.0 X 10"^ < D/H < 3.6 X 10'^, (4.5) 
0.228 <Yp< 0.248. (4.6) 

More specifically, wc impose ^j^^^ < E^sYx where EyigYx is defined in [11]. The largest 
allowed values of E^sYx, as a function of tx, for the assumed abundances oi D/H (low) + 
Yp{IT) can be read out of fig. 43 of [11], for which Bhad = ^-.f^x = 1 TeV and the effect of 
photo-dissociation due to electromagnetic showers was turned off. 

Note that in [11] much less conservative error bars on D/H were used than in [10]. 
Our constraints from HAD showers are therefore going to be accordingly somewhat more 
uncertain than from EM ones, despite not applying the constraint from ^Li at all. 

As we will see, the ranges of supersymmetric parameters excluded by the above con- 
straints will strongly depend on the assumed abundances of the light elements. Below we 
will illustrate this point explicitely in the case of the EM showers by considering the case 

< EyisYx where EyigYx in this case is read off from fig. 42 of [11] where Bhad = 0. 

It would be helpful to have available in the literature bounds on the hadronic fluxes 
assuming more generous inputs, like (4.2)-(4.4), and to allow B^ad to vary. On the other 
hand, the dependence in figs. 38-43 of [11] on mx is rather weak. 

• CMB As originally pointed out in [5] and recently re-emphasized in [27, 28], late 
injection of electromagnetic energy may distort the frequency dependence of the CMB 
spectrum from its observed blackbody shape. At late times of interest, energetic photons 
from NLSP decays lose energy through such processes as 7e~ — > 7e~ but photon number 
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remains conserved, since other processes, like double Compton scatterings and thermal 
bremsstrahlung, become inefficient. As a result, the spectrum follows the Bose-Einstein 
distribution function ^ 

/^(^) = ,i./(.TH._i - (4-7) 

where n here denotes the chemical potential. The current bound is [55] 

< 9 X 10"^ (4.8) 

At decay lifetimes of rx < 4 x 10^^ Q},h'^ sec ~ 8.8 x 10^ sec, where Cli,hP is the abundance 
of baryons, the bound (4.8) leads to a constraint on the EM energy release defined 
previously via the relation 



H = 8.01 x 10^ ^ ' e'^^"^^^^)^'^ X ( ••";\7 ) , (4.9) 



1 sec 



IGeV 



where 



Tdc = 1-46 X 10S(To/2.7K)-i2/5(j^^^2-)4/5(^ _ yp/2)4/5 (4^0) 
~ 6.085 X lOSec (4.11) 
and we have taken Tq = 2.725 if, J^^/i^ = 0.022 and Yp = 0.24. This leads to 

1 /9 

< 1.59 X 10-8 e^r,c/rx)'^^ (|l|C^ . (4.12) 

At later decay lifetimes {tx > 4 x 10^^ ilj/i^ sec ~ 8.8 x 10^ sec), elastic Compton 
scatterings are not efficient enough to lead to the Bose-Einstein spectrum [56]. Instead, 
the CMB spectrum can be described by the Compton y parameter, 4y = Se/e given by 

where T{t) is the CMB temperature and t^s = [r(l — P)]^^^tx, with the Gamma function 
r for a time-temperature relation T oc t~^. In the relativistic energy dominated era in the 
early Universe, for T < 0.1 MeV, 

/ f \-V2 

T = 1.15 X 10"^ GeV, (4.14) 



. 1 sec ^ 

which gives (3 = 1/2. Thus teff = [r{l/2)]'^Tx = tttx- 
The present upper limit on y parameter is [57] 

|y| < 1.2 X 10-^ (4.15) 

This translates into an upper bound on 

< 7.84 X 10-9 (^) ^ 4.42 x lO-'^-L. (4.16) 



Thus we can see that, at the late times Tx ^ 8.8 x 10^ sec, as specified above, the constraint 
on the parameter space coming from the y parameter is applicable while at earlier times 
the //-parameter constraint is applicable [56]. 
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• Collider and Non Collider Bounds The relevant bounds from LEP are from chargino 
and Higgs searches [58] 

m ± > 104 GeV, (4.17) 
nth > 114.4 GeV. (4.18) 

In addition, a good agreement of the measured BR{B — > X57) with a Standard Model 
prediction places strong constraints on potential SUSY contributions to the process, which 
at large Higgs VEV ratio tan (3 can be substantial. We impose [59] 

BR{B Xs-i) = (3.34 ± 0.68) x 10^1 (4.19) 

Finally, we exclude cases leading to tachyonic solutions and those for which the gravitino 
is not the LSP. 

5. Results 

Mass spectra of the CMSSM arc determined in terms of the usual five free parameters: 
the previously mentioned tan/3, mi 12 and mo, as well as the trilinear soft scalar coupling 
^0 and sgn(/i) the sign of the supcrsymmctric Higgs/higgsino mass parameter ^. For a 
fixed value of tan /3, physical masses and couplings are obtained by running various mass 
parameters, along with the gauge and Yukawa couplings, from their common values at 
Mgut down to mz by using the renormalization group equations. We compute the mass 
spectra with version 2.2 of the package SUSPECT [60]. 

We present our results in the usual (m.1/21^0) plane for two representative choices of 
tan/9 = 10 and 50 and for = and > 0. In light of the recent determinations from 
the Tevatron, we fix the top mass at mt = 178 GeV [61]. In fig. 1 we consider the case 
rn^ = 0.2mi/2 while in fig. 2 we take = 0.2mo (top row) and mg = mo (bottom row). 
In both figures we fix the reheating temperature at Tr = 10^ GeV. 

To help understanding the figures, we remind the reader of some basic mass relations. 
The mass of the gluino is roughly given by rrig ~ 2.7mi/2- The mass of the lightest 
neutralino, which in the CMSSM is almost a pure bino, is ~ 0.4mi/2- The lighter stau 
fi is dominated by tr and well above mz its mass is (neglecting Yukawa contributions at 
large tan/3) roughly given by m^^ 2± mp + 0.15m^^2- This is why at mo ^ 'mij^, the stau 
is lighter than the neutralino while in the other case the opposite is true. The boundary 
between the two NLSP regions is marked in all the figures with a roughly diagonal dotted 
line. (In the standard scenario the region of a stable, electrically charged stau relic is 
thought to be ruled out on astrophysics grounds.) Regions corresponding to the lightest 
chargino and Higgs masses below their LEP limits, (4.17) and (4.18), respectively, are 
excluded. Separately marked for tan /3 = 50 is the region inconsistent with the measured 
branching ratio BR{B — > X^^) (4.19). (For tan/3 = 10, and generally not too large tan/3, 
this constraint is much weaker and "hidden" underneath the above LEP bounds.) Here it 
is worth remembering that the constraint is derived by assuming minimal flavor violation 
in the squark sector - the scenario where the mass mixing in the down-type squark sector 
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Figure 1: The plane ('7ii/2, Wo) for tan/3 = 10 (left window) and tan /? = 50 (right window) and for 
Aq = and /j, > 0. The light brown regions labelled "LEP x^" and "LEP Higgs" are excluded by 
unsuccessful chargino and Higgs searches at LEP, respectively. In the right window the darker brown 
regions labelled "6 S7" and the dark grey region labelled "TACHYONIC" are also excluded. In 
the dark green band labelled "figft,^" the total relic abundance of the gravitino from both thermal 
and non-thermal production is in the favored range, while in the light green regions which are 
denoted "NTP" the same is the case for the relic abundance from NTP processes alone. Regions 
excluded by applying conservative bounds on electromagnetic showers from D/H+Yp + '^Li obtained 
with inputs (4.2)-(4.4) are denoted in orange and labelled "EM". Regions excluded by imposing 
less conservative bounds on hadronic showers from D/H + Yp derived assuming Bhad = 1 and 
input (4.5)-(4.6) are denoted in blue and labelled "HAD". (The overlapping EM/HAD-excluded 
regions appear as light violet.) A solid magenta curve labelled "CMB" delineates the region (on the 
side of the label) inconsistent with the CMB spectrum. The cosmologically favored (green) regions 
are ruled out when we apply bounds from D/H + Yp derived with (4.5)-(4.6) as input, as described 
in the text. 



is the same as in the corresponding quark sector. However, even small perturbation of the 
assumption may lead to a significant relaxation (or strengthening) of the constraint from 
6 — > 57 [59]. At small uiq and large tan/3 some sfermions become tachyonic. Finally, for 
some combinations of parameters the gravitino is not the LSP. We exclude such cases in 
this analysis. 

Let us first concentrate on the regions where the total gravitino relic abundance Q^h? 
is consistent with the preferred range (4.1). In all the windows these are represented by 
green bands and labelled (On the left side of the green bands VL^h^ < 0.094 while 

on the other side VL^h? > 0.129.) Their shape looks rather different in both the neutralino 
and the stau NLSP regions. In the former, Q^^h? is mostly determined by neutralino 
decays (NTP mechanism), except that it is relaxed relative to the case of neutralino LSP 
by the mass ratio m^/m^, as in (1.2). (Compare with the standard neutralino LSP case 
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Figure 2: The same as in fig. 1 but for — 0.2mo (top row) and — ttiq (bottom row). In 
the light grey areas the gravitino is not the LSP. Applying bounds from D / H + Yp derived with 
(4.5)-(4.6) as input rules out most of the cosmologically favored regions, except for small patches 
for tan/3 = 50 in the vtLq — 0.2mo and mg windows, as described in the text. 



in, e.g., fig. 5 of [47].) In this region of not too large mi/2 thermal production remains (so 
long as Tr < 10^ GeV) fairly inefficient, since it is proportional to rrig which is still not 
very large there (although growing with mii2). In the right column of the windows where 
tan/3 = 50, one can notice a characteristic resonance due to efficient XX annihilation via 
the pseudoscalar Higgs A exchange. 
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In all the windows light green bands (labelled "NTP") delineate the regions below the 
dotted line (stau NLSP case) for which is consistent with (4.1). These regions 

Or 

become cosmologically favored when one does not include TP or when Tr <^ lOr GeV. 
Again, these are the regions where the relic abundance of the stau (if it had been stable) 
would be consistent with the range (4.1), but shifted to the right and/or up by the factor 
mg/mf^. Note that such regions correspond to ranges of mi/2 beyond those considered 
in [31] and were not identified there. 

However, the contribution from thermal production, which linearly grows with Tr (1.1), 
cannot be neglected when Tr is large. In all the windows of figs. 1 and 2 one can see how 
the effect of TP shifts the cosmologically preferred region from the light green of O^^h? 
to the full green of Q^,h^ where TP dominates. In fig. 1 both rrig and scale up with 
mi/2- As a result, TP dominated regions of O^/i^ are independent of ttiq (compare (1.1))- 
Even though in this figure the contribution from TP is still sub dominant, it does have a 
sizable effect of shifting the vertical bands of total Og/i^ to the left of the ones due to NTP 
production alone. 

In fig. 2, grows with mi/2 (because rUg oc mi/2) but decreases with increasing 

mo (because oc mo)- Compare again (1.1). This causes the green bands to bend 
dramatically towards the diagonal. 

Increasing mg reduces the effect of TP. It simply becomes harder to produce them in 
inelastic scatterings in the plasma. This can be seen in fig. 2 by comparing the bottom 
row (where mg = mo) with the top row (where mg = 0.2mo). The green bands in the 
stau NLSP region, where nig is large, are markedly shifted to the right. (Increasing instead 
mg as a fraction of mi/2 is a less promising way to go as this causes the region where the 
gravitino is not the LSP to rapidly increase.) 

It is thus clear that, so long as Tr < 10^ GeV, one finds sizable regions of rather large 
mi/2 and much smaller mo consistent with the preferred range of CDM abundance. We 
now proceed to discussing constraints from BBN and CMB. 

Constraints from EM showers mainly due to hard photons in neutralino NLSP decays 
into gravitinos (3.2) have traditionally been regarded as severe, and this is confirmed in 
our figures. Even with only the bounds from D/H + Yp + '^Li, derived from conservative 
inputs (4.2)-(4.4), which are labelled as "EM", most of the neutralino NLSP regions are 
ruled out. One exception is when the number density of NLSP neutralinos is reduced, like 
around the "spike" of the A resonance at large tan (3 = 50. On the other hand, in the stau 
NLSP region the constraint from EM showers is in most cases somewhat weaker. 

Adding a constraint from ^Li (not shown in the figures), as in [31], does not actually 
make much difference in all the cases presented. Its main effect appears to be severely 
tempering, to the point of almost removing, the spike regions around the A resonance in 
the neutralino NLSP region. 

On the other hand, adopting the sharper inputs (4.5)-(4.6) into the bounds from only 
D/H + Yp foT constraining electromagnetic showers does lead to a dramatic effect. We will 
discuss this case in more detail below. 

Next we discuss an impact of the constraint from avoiding excessive hadronic fluxes 
(labelled in the figures as "HAD"). Applying the bounds from only D/H + Yp but with 
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Figure 3: Left panel: A comparison of the relative energy {X = x^^i ^-nd * = cm, had) 
with BBN constraints for the case tan/? = 10, /j, > and mg = mo = 62.85 GeV. The black 
long-dashed curves corresponds to and should be compared with the red thick (thin) sohd 
line denoting the upper bounds from D / H + Yp on the relative energy release from electromagnetic 
showers of Cyburt, et ai, [10] (Kawasaki, et ai, [11]), as explained in the text. Clearly, the excluded 
ranges of tx strongly depend on the assumed experimental bounds. The magenta short-dashed 
curves corresponds to C/j^d should be compared with the thick blue solid curve denoting the 
upper bounds from D/H + Yp on the relative energy release from hadronic showers only (no photo- 
dissociation) of [11], as explained in the text. 



less conservative inputs (4.5)-(4.6) (and assuming Bhad = l)i basically rules out the whole 
neutralino NLSP region, thus confirming the conclusions of Feng, et al, [30]. (The cases 
where the hadronic constraint is stronger than the electromagnetic one are marked in blue. 
The opposite case is marked in pink.) It also rules out some cases below the dotted line 
where the very heavy stau NLSP decays fast enough to enhance the importance of bounds 
from hadronic showers. It is, however, possible that, with more conservative inputs, the 
hadronic constraint would not be as severe even in the neutralino NLSP region. 

Last but not least, bounds on allowed distortions of the CMB spectrum prove in 
many cases to be the most severe. They are shown as a solid magenta lines. Regions on 
the side of the label "CMB" are ruled out. While not as competitive in the neutralino 
NLSP region, for stau NLSPs the CMB shape constraint due to the bound on the chemical 
potential (4.8), as already emphasized in [30], and also on the y-parameter (4.16), typically 
provide the tightest constraint. The former bound, being applicable at not very late decay 
times Tx > 8.8 x 10^ sec, excludes regions of the {mi/2,iTT'o) plane closer to the magenta 
curve, while the latter bound, which applies to later decay times, excludes points at smaller 
and/or rriQ. It does not affect the magenta CMB exclusion lines in the ("ii/2)^o) 

plane. 



Coming back to the BBN constraints, in fig. 3 we illustrate a rather strong sensitivity 
of the excluded regions in the (7711/2,1^0) plane to the assumed abundances of the light 
elements. We take tan/3 = 10, = mo and fix mo = 62.85 GeV. In the left panel we plot 
the electromagnetic energy release (black long-dashed line) and compare it with the 
upper bounds taken from fig. 7 of [10] (solid thick red) and with the one corresponding to 
D/H m. fig. 42 of [11] (solid thin red), as described in the previous section. (Constraints 
from "^Li are in this case weaker and are not marked.) Clearly, while the x^NLSP region 
is excluded by both upper bounds, much larger regions of tx are excluded in the fi- 
NLSP region by the solid thin red line. We can clearly see that the excluded regions 
of the (mi/2,'7io) plane strongly depend on the assumed experimental bounds. When 
we apply the EM constraints from [11] to the cases presented here, in most of them the 
cosmologically favored regions (both green bands of and light green ones of il~"'^^/i^) 

become excluded. 

We also plot the hadronic energy release (magenta short-dashed line) and compare 
it with the upper bounds taken from fig. 43 of [11] (solid thick blue line), as described in 
the previous section. The right panel shows tx as a function of mi/2 in order to help relate 
fig. 3 to the lower left panel of fig. 2. 

Given a significant squeeze on the gravitino CDM scenario imposed by the interplay 
of the BBN and CMB constraints, a question arises whether one can find allowed cases at 
Tr even higher than 10^ GeV presented in figs. 1 and 2. As Tr grows, the contribution 
from TP also grows and the green band of the favored range of ri^/i^ moves left, towards 
excluded regions. In all the cases presented in figs. 1 and 2, except for two, even a modest 
increase in Tr, is not allowed by our conservative bounds from BBN and CMB. The two 
surviving cases are presented in fig. 4 for Tr = 5 x 10^ GeV. They are already inconsistent 
with bounds on electromagnetic cascades from D/H + Yp alone when one adopts the less 
conservative inputs (4.5)-(4.6). The case in the right window also becomes excluded by the 
bounds from D/H + Yp + '^Li + ^Li, used in [31], derived using conservative inputs (4.2)- 
(4.4). Finally, an improvement of some order of magnitude in the upper bound (4.8) on 
would also rule these cases out. Given the above discussion, it is unlikely that at higher 
Tr any cases of the favored range of fig/i^ will remain consistent with CMB and/or BBN 
constraints. 

It is interesting that, in all the cases presented above for which mg oc mo, the (pale 
green) region of Jl-^^/i^ ~ 0.1 is almost completely excluded by a combination of BBN 

(-r 

and CMB constraints. In other words, for this case, a contribution to Q,Qh from TP has 
to be substantial in order to escape the above constraints. On the other hand, this is not 
appear to be the case for oc mi/2- 

6. Implications for Leptogenesis and SUSY Searches at the LHC 

It is interesting that, in spite of tight and improving constraints from CMB and BBN, the 
possibility that in models of low energy SUSY with gravity-mediated SUSY breaking the 
gravitino may be the main component of the cold dark matter in the Universe remains 
open. Clearly this is the case when the reheating temperature Tr < 10^ GeV for which a 
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Figure 4: The same as in fig. 1 but for Tr = 5 x 10^ GeV and tan/? = 10 and m^. = mo (left 
window) and tan/3 = 10 and — 0.2mi/2 (right window). Applying bounds from D/H + Yp 
derived with (4.5)-(4.6) as input rules out the cosmologically favored regions, as described in the 
text. 



contribution from thermal production can be neglected. Then the cosmologically favored 
regions due to NLSP decays (non-thermal production) alone shift to larger mij2, which is 
typically (albeit not always) somewhat less affected by the above constraints. 

Perhaps even more interesting is the fact that the reheating temperatures as large as 
5 X 10^ GeV seems to be allowed. This should be encouraging to those who favor thermal 
leptogenesis as a way of producing the baryon-antibaryon asymmetry in the Universe. 

We stress that the above conclusions depend rather sensitively on the assumed input 
from BBN bounds. Given a number of outstanding discrepancies in determinations of the 
abundances of light elements, especially of ^Li and ^He, in our analysis we have decided to 
apply rather conservative bounds derived using rather generous inputs, but also discussed 
the (severe) impact of sharpening them. In constraining electromagnetic fluxes we included 
bounds from D/H + Yp + '^Li derived with conservative inputs (4.2)-(4.4). We did not take 
into account bounds from ^Li nor "^He which otherwise would be most constraining. The 
constraints on hadronic fluxes that we have adopted are somewhat less conservative since 
they are based on less generous inputs (4.5)-(4.6) and on assuming Bhad = 1- Improvements 
in determining the above ranges is likely to provide a very stringent constraint on allowing 
high reheating temperatures Tr ~ 10^ GeV, and perhaps even on the whole hypothesis of 
gravitino CDM in the CMSSM and similar models. Likewise, the (independent from BBN) 
bound from the CMB spectrum, if improved by at least one order of magnitude, is likely to 
prove highly constraining to the scenario. Certainly it would rule out the presented above 
cases of Tr = 5 X 10^ GeV. Finally, we have neglected other possible, although strongly 
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model dependent, mechanisms of producing gravitinos, e.g. from inflaton decays, which 
could also contribute to 

Cr 

An experimental verification of the scenario will, for the most part, be rather chal- 
lenging but not impossible. At the LHC one will have a good chance of exploring some 
interesting ranges of mass parameters. The most promising signature will be a detection of 
a massive, stable and electrically charged (stau) particle. In most cases, the cosmologically 
favored regions correspond to mo < 1 TeV and large ^ 1 TeV, up to 2TeV at small 

tan P, and 4 TeV at large tan (3. A significant fraction of this region should be explored in 
the stau mode since ~ niQ + O.lbmy^ and stable stau mass will probably be probed up 
to ~ 1 TeV. Gluino search may be less promising since mi/2 ^ 1 TeV implies nig > 2.7 TeV, 
which will be just outside of the reach of the LHC. However, some interesting cases may 
still be explored. One is that of fairly small tan/3 and oc ttiq, as in the left windows 
of figs. 2 and 4. The other is that of large tan f3 and large Tr, as in the right window of 
fig. 4. On the other hand, since in the cosmologically favored regions tuq <C mi/2 squarks 
(and sleptons) will be fairly light and lighter than the gluino. More work will be needed 
to more fully assess to what extend the cosmologically favored regions will be explored at 
the LHC. Finally, the staus will eventually decay. The proposal of [62] of measuring their 
delayed decays is in this context worth pursuing as it would give a unique opportunity of 
experimentally exploring the hypothesis of gravitino cold dark matter and of probing the 
Planck scale at the LHC. 
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